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Abstract 

In this paper, we propose observers for the dynamics of anthracnose disease described in |12] . Spatial 
and non spatial versions of the observers are given following an approach similar to |12] . The models 
given in m are improved in order to be more realistic. The conditions on parameters are more general. 
There are also changes in the equations modelling the dynamics of the berry volume (v) and the rot 
volume (vr). We study theoretically the proposed observers in terms of well-posedness and convergence. 
We make several simulations to assess the effectiveness of observers which definitely display fairly good 
behaviour. 
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1 Introduction 

Anthracnose is a phytopathology which occurs on several commercial tropical crops. The coffee is 
concerned by that disease due to the Colletotrichum kahawae which is an ascomycete fungus [SlElEllia 
Several models of the athracnose dynamics have been proposed in the literature p nni ttn ng 
[izi uBi ng Eg [21] in order to better understand and control it. Recently, in [ig[i 3 ] the authors studied 
the optimal control of a diffusion model of anthracnose with continuous and impulsive strategies. Optimal 
controls were computed with respect to given cost functionals. The general model proposed in [T2| was given 
by the following equations 


dtd = a (t, x) (1 — re (t, x) 6) + div {A (t, x) VO) , on x U (1) 

{A {t, x) VO {t, x ), n (x)) = 0, on x dU (2) 

where n (x) denotes the normal vector on the boundary at x and 

0 (0, x) = p (x) € [0,1[, X G U C (3) 

dtV = j3 (t, X, 0) (1 - v/ ((1 -0)7] {t, x) Umax)) (4) 

dtVr = ^{t,X,0){l-Vr/v) ( 5 ) 

{v (0, x) , Vr (0, x)) € ]0, Umax] X [0, Umax] , X G U C (6) 

where tc (t, x) = 1 / (1 — au {t, x)) and v (0) > Vr (0). 


In the model above, 0 is the inhibition rate, v and Vr are respectively the fruit volume density and the 
rot volume density. The function v is upper bounded by a value Umax which models the natural fact that 
the fruit growth is limited. Nonnegative functions a, /3 ,7 : R+ x [/ x R ^ R_i_ characterize the effects of 
environmental and climatic conditions on the rate of change of inhibition rate, fruit volume, and infected 
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fruit volume respectively piiniiii]. There is a control u representing the chemical strategy consisting 
on the effects after application of fungicides. The positive term 1 — a € ]0, 1[ models the inhibition rate 
corresponding to epidermis penetration by hyphae. Once the epidermis has been penetrated, the inhibition 
rate cannot fall below 1 — cr, even under the maximum control effort (u = 1). In the absence of control effort 
(n = 0), the inhibition rate should increase towards 1. The effects of environmental and climatic conditions 
on the maximum fruit volume are represented by the parameter r/ which is a ]0, 1] —valued function. There is 
a spatial spread of the disease in the open domain [/ C of class through the diffusion term div {AV0). 
The boundary condition {A'V9,n) = 0 where is a 3 x 3-matrix (aij) could be viewed as the law steering 
migration of the disease between U and its exterior. For instance, if A is reduced to the identity matrix I 
then (V9,n) = 0 means that the domain U does not have any exchange with its exterior. 

Also in m a within-host model has been studied and given by the following equations : 


dt9 = a{t) {1 — w {t) 9) (7) 

dtv = 13 {t, 9) (1 - v/r] (t) Umax (1 - 9)) (8) 

dtVr = -f{t,9){l-Vr/v) (9) 

{9 (0) ,U (0) ,u,. (0)) G [0, 1[ X ]0,Umax] X [0, Umax] (10) 


where w (t) = 1 / {1 — au (t)) and u ( 0 ) > Vr ( 0 ). The parameters in ([7]) — ffTOjl still have the same interpre¬ 
tations given in the general model ([I]) — dSD- 

In several cases, especially for the results in na m on anthracnose disease, an optimal control 
strategy is given as a feedback and needs to know the state of the system and parameters values. It could 
be difficult to know exactly the dynamics of the inhibition rate while it is easier to observe volumes u 
and Vr- The aim of this paper is to provide some observers for the state 9 in the models ([I]) — dS]) and 
(0 ~ ijlOP assuming that volumes u and Vr are known. An observei0 is a function 9 having its own dynamics 
depending on the available observations (u and Vr) and converging towards a neighborhood of 9. Let us 
recall the following useful adapted definitions about the stability of an observer. The authors refers to 
[El na Ea |23] and references therein for more general and classical definitions. 

Definition 1.1 Assume that 9 is S-valued. An observer 9 is said to be 


(i) (Locally) stable if there are two neighborhoods U,V C S ofO and a time to > 0 such that if 9 {0) — 9 (0) € 
U then Vt > to, 9 (t) — 9 (t) € V; 

(ii) (Locally) asymptotically stable if there is a neighborhood U S of 0 such that if 9 (0) — 0 (0) € 1/ then 
lim 9 {t) — 9 (t) = 0; 

t^OO 


{Hi) (Locally) exponentially asymptotically stable if there are a neighborhood U C S of 0, a positive real 
constant k and a time to > 0 such Vt > to, 9 {f) — 9 (t) < exp {—kt). 


Those local stabilitiy properties are said to be global if one can extend the set U to the whole set S. 


The remainder of the work is organized as it follows. In section [2] we present some new realistic 
considerations made for the models. We also show that the new model is well-posed. We design observers 
for the non spatial model in section [31 A theoretical study is done in subsection 13.11 and we make numerical 
evaluations in subsection 13.21 We also design observers for the spatial model in section [3] and we study them 
theoretically in subsection 14.11 A numerical evaluation is carried out in subsection 14.21 Finally, there is a 
conclusion in section [5l 

^See for more details on the topic. 
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2 New modelling of anthracnose dynamics 


In this section we introduce some changes in the models given in [12] based on some remarks in order 
to become more realistic. Our modifications are especially made on the description of the dynamics of v 
and Vr- 

Looking at equation (|3|) we can see that when 0 tends to the value 1, v tends with an infinite speed to 
zero. However, we think that either there is a disease or not the berries will reach and remain greater than 
a minimal volume. Indeed, the disease can not be triggered when the berry is less than a minimal size as 
for instance the bud size. Even if the fruit totally rots it keeps a minimal size. We then suggest to add a 
very small nonegative term e <C 1 in such way that equation Q becomes 

dtV = 13 {t, X, 0) (1 - v/ {{1 + e - 9)r] {t, x) Umax)) (11) 

With the same arguments equation dHj) becomes 

dtv = 13 (t, 6) (l-v/r] (t) Umax (1 + e - 9)) ( 12 ) 

The use of the term e generalizes the model in such way that we can recover the original model given in 
m by setting e = 0. With a positive e even if 0 = 1, u might be also positive. Indeed, the berries reach a 
minimal volume in the neighborhood of eumax™iii{?/ (t)}. 

We now focus on the equations dSj) and dHj). In [T2| it was established that Vr is [0, Umax] valued. However, 
that behaviour is not sufficient. Indeed, Vr should remain less than u. On the other hand Vr should remain 
equal to zero while u is less than the minimal value corresponding to the threshold volume where the disease 
can be triggered. We then suggest that Vr remains equal to zero when u is not greater than 0. Instead 
of studiying directly u^, we introduce a proportion function p depending on time and space variables, and 
given such as 

Vr = pv (13) 

We propose the following dynamics for the rot proportion p: Vx € U, p (0, x) G [0,1] and 

dtp (t, x) = 7 (t, X, 9 {t, x) ,v{t,x) ,p {t, x)) (1 - p (t, x )), V {t, x) € x U, (14) 

where 7 ; M+ x [/ x ^ R is a function satisfiying V (t, x) G R_|_ x U, 

Assumption 2.1 j is a measurable with respect to the two first parameters and locally Lipschitz continuous 
with respect to the third parameter. 

Assumption 2.2 7 (f,x,.,.,0) is nonnegative and such that 

7 (t, X, 0,., 0) = 0 = 7 (t, X, ., 0,.) . (15) 

Assumption 2.3 7 (t, x,.,.,.) is increasing with respect to the first parameter and 7 (t, x, 0,.,.) is decreasing 
with respect to the last parameter. 

For instance 7 could be chosen with the simple but general form 7 (t, x, yi,y 2 , 2 / 3 ) = (71 {t, x) yi — 72 {t, x) 2 / 3 ) y 2 , 
with 7 i, 72 : R+ X [/ X R^ —>■ R_|_. If we omit the dependence with the space variable we get the corresponding 
equation for the non spatial model: p (0) G [0,1] and 

dtp (t) = l{t,9 (t) ,v{t),p (t)) (1 - p (t)), Vt G R+. (16) 

As we will check in the sequel, the function p is [0, l]-valued. The assumption (12.211 guarantees that 
while the berry has a null volume (without berry) or the disease has not started the rot volume remains 
null. The assumption (12.311 means that the rot volume increases with the inhibition rate. When there is a 
not inhibition the volume of rot does not increase while the fruit grows better and therefore the proportion 
p decreases. 
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2.1 Well-posedness of the problem for the non spatial model 

This subsection is devoted to the proof of the well-posedness of the non spatial model. We make the 
following useful assumptions. 

Assumption 2.4 a € (M+; [0,1[). 

Assumption 2.5 rt, ry G L°° (M+; [0,1/ (1 -|- e)]) and Vt > 0, ry (t) > ry* G ]0,1[. 

Assumption 2.6 13 : M+ x [/ x M —M_|_ is nonincreasing function with respect to the third parameter. 
Vz G M, /3 (., .,z) is a measurable function and belongs to the space (R+ x R;M). 

Assumption 2.7 Vt > 0, the function f3 ft, .) is differentiable and dxj3 {t, .) G (R;R). 

The assumptions ()2.4ll - (l2.6p are more general than those given in [12] in the sense that all the parameters 
of the model were assumed continuous with respect to the state variable while here they are just measurable 
and essentially bounded. The following existence and uniqueness propositiont holds: 

Proposition 2.1 If {6 (0), r (0), p (0)) G [0,1] x [0, Umax] x [0,1] then the problem d?]), ffT^ and (fTHIl has a 
unique absolute continuous solution {9,v,p) valued in the set [0,1] x [0, Umax] x [0,1]. 

Proof. Using the Caratheodory theorem (see m) we get easily the existence of an absolutely continuous 
solution. The logistic form of the equations d?]) and (I12p shows clearly that for a maximal solution {9,v,Vr) 
the couple (0,u) is necessarily valued in the set [0,1] x ]0,Umax]- On the other hand from the equation (fTHD . 
p is valued in [0,1]. To show it let set pi = max{0, —p} and p 2 = max{0,p — 1}. The functions pi and p 2 
are continuous as p is continuous. Without loss of generality we will assume that pi and p 2 are positive on 
an open interval of time ]0,T[. Then Vt G ]0,T[, 

Pi{t) = - [ 7 (s, 9{s),v (s), -pi (s)) (1 Pi (s)) ds (17) 

Jo 

< 0 

P 2 {t) = -[ j{s,9{s),v{s),p2{s) + l)p 2 {s)ds (18) 

Jo 

< 0 

There is a contradiction meaning that pi and p 2 are identically null. Therefore, p is [0, l]-valued. ■ 

2.2 Well-posedness of the problem for the spatial model 

In this subsection the well-posedness of the spatial model is proved. To that aim, we make some 
assumptions. 

Assumption 2.8 a G L“^(R+;L°° ([/;R+)). 

Assumption 2.9 u,ri € (R+; (U; [0,1])) and Vt > 0, 1/ (1 -|- e) > ry ft,.) > rj* ^ L°° {U; [0,1]) with 

0 < inf (ry*) 

Assumption 2.10 (3 : R+ x [7 x R —>■ R+ is nonincreasing function with respect to the third parameter. 
Vz G R, /3 (., .,z) is a measurable function and belongs to the space (R+ x R;R+). 

Assumption 2.11 Vz,j G {1,2,3}, Oij G (R+; (t/; R)). 
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Assumption 2.12 3(^ € such that Vt G M+, Vw € ([/;] 


/ {A{t,x)\7w{x),Vw{x))dx>C (\7w (x) ,\7w (x)) dx. 

Ju Ju 

Assumption 2.13 Vt > 0,Vx G C/, i/ie function /3 (t, x,.) is differentiable and dyj5 {t,x ,.) G 


Assumptions ()2.8p - fl2.10p are also more general than those given in |12| in the sense that all the parameters 
of the model were assumed continuous with respect to the state variable while here they are just measurable 
and essentially bounded. From assumptions (I2.1ip - (l2.12p the following problem has a unique solution in 
(U) for an arbitrary but fixed time t > 0. 

div (A (t, x) S/w (x)) = / (x), Vx G [/ 

(A {t, x) Vrc (x), n (x)) = 0, Vx G dU 

where / G {U). Following [3] in Theorems 3.6.1 and 3.6.2, there is an orthonormal complete system 
{ipn {t, C (U) of eingenfunctions and eingenvalues {A„ (t)} such that Vn G N, 

div (A (t, x) V(pn (t, x)) = An (t) (Pn (t, x), Vx G f7 
(A (t, x) V(pn (t, x) , n (x)) =0, Vx G dU 

Moreover, {ipn {t, ■)}ngN (^) class then ((/?„ (T ■))nGN (^) —valued. 

Now we make the following additional assumption. 

Assumption 2.14 The sequence {<Pn) does not depends on the time (ie ipn {t, .) = (fn (■), Vt > Oj. 


Assumption (12.1411 could happen if A (t,.) has the form p, (t) B (.) with p, (t) G M, Vt > 0. It will be the 
case in particular if A (t, ■) = p (t) I and therefore div (A [t, .) S/w) = p, {t) Aw. Here I denotes the identity 
matrix of M^. Whether (12.1411 is satisfied a weak solution 0 of ([T]) — (l3|) can be written as the sum > OnTn 

• ^71=0 

where each 9n is an absolutely continuous function of the time and satisfies Vt > 0, 

E OO P 

I Tm (x) qPn (x) dx (f9) 

m=0 Jjj 

= (0) + / Xn (s) On (s) ds + / / a {s, x) ifn (x) d.sdx 

Jo Jo Ju 

— a {s,x) w {s,x) 6 {s,x) ifn (x) dxds. 

Jo Ju 

Proposition 2.2 The problem ([I]) — ([3]) has a unique solution 6 ^ C (M+; ([/; [0,1])) which is absolutely 

continuous with respect to the time. 

Proof. It is sufficient here to etablish that if there is a maximal solution then that solution is valued in the set 
[0,1]. Indeed, conditions (|2.8p - (l2.12ll are sufficient to use the Caratheodory theorem (see in [T] the methods 
used in the proof of the Lemma A 2.7, pl87-191). Let Vt > 0, = max{0, —9} and 02 = max {0,0 — 1}. 

We have 0i (0,.) = 02 (0,.) = 0. Let Si be a convex open subset of M+ where functions 0* (t,.) (i G (1, 2}) 
are positive on subsets of U with a positive measure. Then Vt G Si we have 


dt ll/i (L •)IIl2(c/;R) 


u 

= - a{t,x){l + w {t, x) fi {t, x)) fi {t, x) dx 

Ju 

< 0 
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||/2(t, .)llL 2 (r/;K )/2 = / f 2 {t,x)dtf 2 {t,x)dx 


u 


/ a {t,x) {1 — w {t,x)) f 2 {t,x) dx — / a {t,x) w {t,x) f 2 {t,x) dx 

Ju Ju 


< 0 

Using the above inequalities we get that functions 9i {t, .) are nonpositive on the sets Si which are necessary 
empty. We conclude that Vt > 0, 6 {t, .) is valued in [0,1]. ■ 

Proposition 2.3 The problem (ITT]) . (|14|) has a unique solution {v,p) G C (M+; ([/; [0, Umax] x [0,1])) 

which is absolutely continuous with respect to the time. 

Proof. Since 9 is valued in the set [0,1] and there is not diffusion in the equations (llip and (11411 . we can 
fix the space variable and use a proof similar to the one given for proposition 12.11 ■ 


3 Observation for the within-host model 

3.1 Theoretical design of observers 

The main objective of this subsection is to design observers for the problem ([7]), (I12p and (|16ll in order 
to estimate the inhibition rate 9. Let us consider the following system : Vt > 0, 


dt9 (t) = a (t) [I — w (t) 9{t)) + ki (t) 4>i [t,9 (t) ,v{t)) + k 2 (t) 4>2 {t,9 {t 


and 


where 


(f>i it,x,y) = 

4>2 {t,x) = 


dtv {t) = 13 \ t, 9 {t)j (j )3 [t, 9{t) ,v (t) 

{l-v (t) /y) {l + e-x),\iv{t)<y and (x, y) € ]0,1[ x 


( 20 ) 

( 21 ) 


and 


0, otherwise 

dtp {t) - 7 (L X, v{t),p (t)) {1- p (t)), if X € ]0,1[ 
0, otherwise 

03 {t,x,y) = 1 - 


{l + e-x)r]{t)vn 


We make the following additional hypothesis : 


Assumption 3.1 ki,k 2 € L 




By a solution of the system (|20ll — (|2ip we mean an absolutely continuous function of the time {9,v 
which satisfies Vt > 0, 

9 {t) = 9 {0) + J a (s) — w (s) 9 (s)^ + J ki (s) cpi [s,9 (s) ,v (s )) ds 

+ J k2{s)4>2 {s,9{s)^ds 


and 


v{t) = v{0) + J P (s, 9 (s)) (/)3 {s,9 (s) ,v (s )) ds. 
We have the existence and uniqueness 
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Proposition 3.1 If ^0(O),i;(O)j € [0,1] x [0.t> mav ] then the problem ([20]) — ([2T]) has a unique solution 
valued in [ 0 , 1 ] x [ 0 , Umax] • 


Proof. We first show that if is a local solution of the system (EOD - dUD then [^iVJ is valued 

in [0,1] X [0,Umax]- Let Vt > 0, 0i = max|o,— 0 |, 62 = max |o, 0 — l| ,ui = max{0,—u} and V2 = 
max{0,u — Umax}- We have 6i (0) = Uj (0) = 0, Vf € {1,2}. Without loss of generality we can assume that 
the functions 9i and Uj (i G {1,2}) are positive on an open interval ]0,T[. Then we have 

9i{t) = — [ a (s) (1 + w (s) 9i (s)) ds — f ki (s) (j)! {s, —9i (s) , u (s)) ds 
Jo Jo 

- k 2 (s) (j )2 {s, -9i (s)) ds 

Jo 

< 0 


and 


a (s) {I -w{s) -w (s) 92 (s)) ds + f ki (s) 4>i {s, 1 + 6*2 (s), u (s)) ds 

Jo 

+ / A:2 (s) (^2 (-s, 1 + ^*1 (-s)) d-s 

Jo 

< 0 . 



That is 9i and 02 are identically null and 9 is [0, 1]-valued. In the same manner one can show that u is 
valued in [0,Umax]- 

Now let show the existence and the uniqueness of the solution of (PUI) — (piD . It suffices to establish 
existence and uniqueness of a local solution and use the fact every local solution of ([2n]l — (1211) is bounded 
to conclude using Theorem 5.7 in [1]. 

Let consider the function F : {t, x, y) G M+ x [0,1] x [0, Umax] 1®^ defined by 

(t, X, y,z) = a {t) (1 - xw (t)) + kicfi (t, x, y) + k 24>2 (L x ), 

and 

(L x,y) = /3 {t, x) (fs {t, X, y) . 

The function F is measurable with respect to the time t and continuous with respect to {x,y). Using the 
Caratheodory theorem there is a local solution of (|20p — ([211) . Moreover, F is Lipschitz continuous with 
respect to (x, y, z). Therefore, the solution is unique and global. Note that the hypothesis (EZD implies that 
Vt > 0, the function /3 (t, .) is Lipschitz continuous and bounded on the set [0,1]. ■ 

Now we state the main results of this subsection. Let define the function 5 : M —{0,1} such as 

,5 (x) = / ® ^ 

^ ' { 0 , otherwise 

Theorem 3.2 Let consider the system (1201) — (12ip and assume that u (0) = u (0). 

(i) If the functions ki and k 2 are identically null and 

inf {a(t) > 0} > 0 (22) 

then 9 is a globally exponentially asymptotically stable observer for 9. 
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(a) If the function ki is identically null and k 2 is not identically null then 9 is a a locally asymptotically 
stable observers for 9. Moreover, if there is a positive function (M+;M+) such that 


It (i, 2/1,2/2,2/3) - 7 (^,^;i,^2 ,^;3)| > (t) ||y - z\ 

then 9 is a a globally exponentially asymptotically stable observers for 9. 
(Hi) If the function k 2 is identically null and Vi > 0, 


(23) 


aw + ki6 ( 0 ) + 


ki6 [ 9 ) f3 {t,9) (ryt^max {l + s-9)-v) 


> 0 . 


aVVVma.^ (1 - Ow) 

then 9 is a a locally exponentially asymptotically stable observers for 9. 

{iv) If the functions ki and k 2 are not identically null, and Vi > 0, 

ki5 (t) 15 ft, 9) fqvrai,^ (1 + e - 0) - i;) 


(24) 


aw + ki5 + k 2 (i) ( 1)2 (t, o'j + 


arjVVma.y, {1 - w9) 


> 0 


(25) 


then 9 is a a locally exponentially asymptotically stable observers for 9. Moreover, if the condition 
is satisfied and 


k 2 (j)2ft,9j > ki(f)i \^t,9,v 
then 9 is a a globally exponentially asymptotically stable observers for 9. 
Proof. Let e = 9 — 9 he the error of estimation. 

(i) The error e satisfies 

dte = —a (i) w (i) e 

and ^ 

e (i) = exp J a{s)w (s) e ( 0 ). 

The result follows. 

{ii) The error e satisfies 

dte^ = —2a (i) w (i) (i) — 2k2 (i) 02 (t,9 (i)') e (i) 


(26) 


(27) 


(28) 


= -a{t)w{f)e^ ft) - 2 k 2 ft) (1 - pft)) ( 7 (i, 0 (i), i; (i), p (i)) - 7 (i, 0 (i), u (i), p (i)) ) eft) 
< —a ft) w ft) e^ ft) 

e^ (i) < exp ( — / a(s)w (s) ds ] ( 0 ) 


and 


Moreover, if the condition (I23p is satisfied then 


and 


< —a ft) w (i) ft) — 2 k 2 (i) (1 — p (i)) ft) e^ (i) 


e^ (i) < exp ( — / {a{s)w (s) + 2 k 2 (s) (1 — p (s)) (s)) ds ) ( 0 ). 






{Hi) We have 


^ 13 {s, 9) _ dtvde/3 {s, 6) - P (s, 9) dedtv 

f3{s,9) — dtv {(3 {s, 9) — dtv)'^ 

_ (3'^ {s,9) {{1 + e - 9) dev + v) 

(1 + e - 9f {13 {s, 9) - dtvf 

and at the neighborhood of 0, 


13 [s,9 


f3{s,9) 


_ 13"^ {s,9) {{1 + e - 9) dev + v) 

(3 (^8,9^ - dtv {Hs,Q)-dtv) r^vra!,^{l + e-9)^{^{s,9)-dtvy 

jyt’max (1 +e - 0) ^■Vmax^ TlVraa.^ {I + E - 9) deV 


-e — 


■e + eO (e) 


e + eO {e) 


It follows that 


V 

V 


/3(s,9 


^■Wmax (1 + £ 9) p q'J _ 


= 1 - 


' +^ie- 


1 + £ — 9 V 


ve 


r/'Vmax {I+E-9) 


0{e). 


finally, 


dte = —a {t) w {t)e — (t, 9, 

ki {t) 6 ( 9 ') (1 + e — 9) dev 


e — [a {t) w {t) — ki {t) 5 [9) )e 


ki {t) 6 (9) V 


Vj^rasux 


-eO (e) 


ki{t)5[9] 13 (t, 9) (r?fmax {l + £-9)-v) 


ariWmax (1 - 9w) 


e — (a{t)w{t) + ki{t)6(9]]e 


The given approximmation on the dynamics of e shows clearly that 0 is a locally exponentially asymp¬ 
totically stable observers for 9 if the condition (I24p is satified. 

{iv) With the same arguments developped in {ii) and {iii) we easily get that 0 is a locally exponentially 
asymptotically stable observers for 9. Using the conditions (I23p and (I26p we are sure that even the 
behaviour of the term with (pi tends to introduce unstability, the stability is maintained by the term 
with (p 2 . 


3.2 Numerical evaluation of the observers 

This section is devoted to numerical simulations in order to check the effectiveness of the observers we studied 
theoretically in the subsection 13.11 The Theorem 13.21 guarantees asymptotical stability of the observers, but 
pratically the convergence is needed in a finite time. Since the coffee cultivation campaign is a cyle of one 
year, our simulations will cover the year. Parameters are taken following |12j . The control strategy u is 
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taken with many variations in order to study their impact on the.performances of the observer. For every 
time t > 0 , tt is given by 

u{t) = sin^ {t — exp (^-uj 2 {t - • (29) 

The functions a, j3 and 7 are taken with the following form 

a (t) = Pi (t) + 61 (1 — cos (cit)) {t — di)^ , Vt G M+, (30) 

/3 (t, x) = 62 (1 - cos (c 2 t)) {t — ^ 2 )^ P2 {x) , V {t, x) € M+ X M, (31) 

and 

7 {t, Xi, X 2 , X 3 ) = 63 (1 - cos (cst)) {t - d'if (xi - KX 3 ) X 2 , V (t, x) G M+ X (32) 

Pi is a nonnegative function of the time t and p 2 is a real positive function of x. Vi G {1,2,3}, bi, Ci, 

and di are positive coefficients corresponding respectively to the maximal amplitude, the pulsation and the 
global maximun of a, /3 and 7 . k is a positive constant regulating the evolution of the rot volume with 
respect to the inhibition rate. The terms 1 — cos{{cit) represent the seasonality probably due to climatic 
and environmental variations. 

Theoretically the functions ki and /c 2 might take abitrary positive numbers and we can set them constant 
without loss of generality. However, for the stability of the numerical scheme we were constrained to set 
them less than 1/lOAt. Here At denotes the constant time step. The initial conditions are taken such as 
0(0) G {0.05,0.25,0.50,0.75}, v (0) G {0,0.25,0.50,0.75}, p(0) G {0.25,0.50,0.75} , 0(0) = 0, v{0) =v{0). 
Indeed, at the beginning the fruits are small and the inhibition rate is relatively small. 

The following table gives the assumed parameters values. 


Parameters 

Values 

Source 

Parameters 

Values 

Source 

bi 

5 In(10) 

Assumed 

Pi (t) 

0 

Assumed 

b2 

Umax In (lO^Umax (1 “ CP*)) /2 

Assumed 

P2 (x) 

2 — X or (2 — x)^ 

Assumed 

63 

Umax In(lO^nmax) 

Assumed 

K 

1 

Assumed 

di 

7.5 X 10-1 

m 

^max 

1 cm^ 

Assumed 

d2 

7.5 X 10-1 

m 

e 

10-^ 

Assumed 

ds 

7.5 X 10-1 

m 

v{t) 

1/ (1 + e) cm^ 

Assumed 

Cl 

IOtt 

m 

a 

0.9 

m 

C2 

lOvr 

m 

h 

0 or 10^ 

Assumed 

C 3 

lOvr 

m 

k2 

0 or 10^ 

Assumed 

Wl 

257r 

Assumed 

At 

10-^ 

Assumed 

0 J 2 

10 

Assumed 

Pi 

0.6 

Assumed 

Pi 

0.4 

Assumed 





Table 1: Simulation parameters for the non spatial model 

In the remainder of the subsection we present simulations of the observers and their estimation relative 
errors dynamics. The relative error seems to be a good mean to evaluate the performance. For each scenario 
we display two groups figures. The first one represents the dynamics both of the inhibition rate and observers 
corresponding to each values of p (0) <0 (0). The second group of figures shows relatives errors of observers 
corresponding to each values of p (0) < 0 (0). 
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3.2.1 Simulations with 0 (0) = u (0) = 0.05 

The figures [T] and [2] show the behaviour of the observers when the observation is started early. As we can 
see the global performance of each observer is relatively good but the case fci = 0 and k 2 = 10^ seems to be 
the better one. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and k 2 = 0 


ki = 10^ and k 2 = 10^ 


Figure 1: Inhibition rate and observers for 0 (0) = u (0) = 0.05. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 ki = 10^ and /c 2 = 10^ 

Figure 2: Relative absolute estimation error for 9 {fi) = v (0) = 0.05. 
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3.2.2 Simulations with 6 (0) = 0.05 and v (0) = 0.5 

The figures El and H] correspond to an observation starting at the beginning of the disease but with a well 
developped berry. We have better performances when k 2 = 10^. With ki = 10^ and A :2 = 0 we notice an 
unstability behaviour in figure 0] although the error takes the null value. That may be due to the constant 
sign of the function 4 >i. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and k 2 = 0 ki = 10^ and k 2 = 10^ 

Figure 3: Inhibition rate and observers for 9 (0) = 0.05 and v (0) = 0.5. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 


ki = 10^ and /c 2 = 10^ 


Figure 4: Relative absolute estimation error for 9 (0) = 0.05 and v (0) = 0.5. 
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3.2.3 Simulations with 6 (0) = 0.75 and v (0) = 0.05 

In the figures [S] and El the inhibition rate is already high while the berry is little. As expected the worst 
estimation corresponds to ki = k 2 = 0. The best estimations are still made when k 2 = 10^. However, when 
ki = 10^, the same unstability behaviour in figure 0] appears. Note that the performances of the observers 
are better with p (0) small. 


= ^2 = 0 


/si = 0 and k 2 = 10^ 


ki = 10^ and A :2 = 0 ki = 10^ and k 2 = 10^ 

Figure 5: Inhibition rate and observers for 9 (0) = 0.75 and v (0) = 0.05. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 


ki = 10^ and /c 2 = 10^ 


Figure 6: Relative absolute estimation error for 6 (0) = 0.75 and v (0) = 0.05. 
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3.2.4 Simulations with 6 (0) = 0.75 and v (0) = 0.5 

For the figures [7] and [5] the inhibition rate is already high and the fruit is relatively mature when the 
observation starts. The worst case are still when k 2 = 0 and the best are with k 2 = 10^. Again the 
performances of the observers are better with p (0) small. It seems that starting observation with 9 (0) small 
and V (0) big guarantees better results. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and k 2 = 0 ki = 10^ and k 2 = 10^ 

Figure 7: Inhibition rate and observers for 6 (0) = 0.75 and v (0) = 0.5. 
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fci = /C2 = 0 


ki = Q and k 2 = 10^ 


ki = 10^ and /c 2 = 0 


ki = 10^ and /c 2 = 10^ 


Figure 8: Relative absolute estimation error for 9 (0) = 0.75 and v (0) = 0.5. 
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4 Observation for the spatial model 


The aim of this section is similar to the previous section. We design observers for the spatial dynamical 
system (HI) - © , © , (HH) and (|14p in order to estimate the inhibition rate 9 and we improve them using 
computer simulations. 


4.1 Theoretical design of the spatial model observers 

Let consider the following systems : Vt > 0,Vx € U, 


where 


4>i {t,x,y,z) 


dtO = a {t,x) — tc {t, x) o'j + Ki (t, x) 4>i x, 9, 

+ K 2 {t, x) 4>2 (t, X, 9 'j + div (^A {t, x) V9^ , on x U 

^4 {t, x) V01 (t, x) , n (x)^ = 0, on x dU 
dtv = P (t, X, 0^ 4>3 (t, X, 9, 

( {1 — V {t, x) jz) {1 + e — y) AA V {t,x) < y and (x, y,z) ^ U x ]0,1[ x 
I 0, otherwise 


f dtp {t, x) - 7 {t, y, V {t, x ), p {t, x)) ( 1-/9 {t, x)), if (x, y) e U x ]0,1[ 
( 0 , otherwise 


and 

$3 {t,x,y,z) = 1 - 
We make the following assumption: 


z 

y) p {t, x) Umax ’ 


(33) 


(34) 

(35) 


Assumption 4.1 Ki,K 2 € (M+; L°° ([/; M+)). 


By a solution of the system pSIl — we mean an absolutely continuous function with respect to the 
time which satisfies Vt > 0,Vx € U, 

9{t,x) = 0(0, x) + I a (s, x) ^1 — w (s, x) 9 (s, x)^ ds + div ^4 (s, x) V0i^ ds 


+ 


[ Ki{s,x)' 
Jo 


4>i s,x,0,u + / A 2 (s, x) 4'2 s, X, 0 


V {t, x) =v 


( 0 ,x) + [ 
Jo 


/3 ( s, X, 0) $3 ( s, X, 0, u ) ds 


We have the following 


Proposition 4.1 If ^0 (0,.), u (0, .)^ G L^(f7;[0,1[ x [0, Umax]) then the system (f33]l — (l35]) has a unique 
weak solution ^0, u^ G C (M+; (U ; [0,1])) x C (M+; (U ; [0, Umax])) • 
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Proof. We first show that if {^9, vj is a local solution of the system (l33|) — (1351) then Vt > 0, \^9 (t, .) , v {t, 

is valued in [0,1] x [0,Umax]- Let Vt > 0, 0i = max|o,—02 = max|o,0 —l|, ui = max{0,—u}, 

V 2 = max {0, V — Umax}- We have 9i (0, -) = 02 (0, -) = ui (0, -) = U 2 (0, -) = 0. Let Sf and be convex open 
subsets of M+ where respectively functions 9i {t, .) and Uj (t, -) (i € {1,2}) are positive on subsets of U with 
a positive measure. Then Vt € Sf (respectively 5}) we have 

2 


dt 


/2= f 9i{t,x)dtfi{t,x)dx 
L^{U-,R) Ju 


/ a{t, x) 

Ju 


1 + w {t, x) 01 (t, x )) 01 (t, x) dx 


lU 


— / Ki {t,x) lt,x,—9i {t,x) ,v] dx 


- J K2{t,x)^2(t,x,-9i{t,x)) dx 
< 0 


O 2 {t,.] 


L2([/;]E 


/2 = / 02 {t, x) dt92 {t, x) dx 


/ a {t,x) {1 — w {t,x)) 92 {t,x) dx — / a {t,x) w {t,x)9^ {t,x) dx 

Ju Ju 


< 0 


dt 11^1 (L ■)lli 2 (f/;K) /2 = {t,x)dtvi {t,x)dx 

= - j l3{t,x, -ui {t, x)) $3 (^t, X, 01 , -ui {t, ui {t, x) dx 


< 0 


and 


dtHh(t,.)Hl2(iy;R)/2= / V 2 {t,x)dtV 2 it,x)dx 


u 


J /? it,X,V 2 {t,x) + Umax) ^3 (t,X,9i,V2 {t,x) + Umax) V 2 (t, x) dx 


< 

< 0 

Using the above inequalities we get that functions 9i {t,.) and Vi {t, .) (i G (1, 2}) are nonpositive on the sets 
Sf and which are necessary empty. We conclude that Vt > 0, ^0 (t, .), u (t,.)) is valued in [0,1] x [0, Umax]- 
Now let show the existence and the uniqueness of the solution of (f33l) — {35]). It suffices to establish 
existence of a local solution to conclude the result. Just as decomposition {T9]l a solution to (f33|) — {Ml) can 

E OO — — 

9n^n- Each 9n satisihes the realtion 

r).=n 


'— ^ ^ oo 

dn (^) — / ndm I (tt') ( 2 ^) dx 

^m=0 Jjj 

= 0n (0) + / Xn{s)9n{s)ds + / (fn (x) a {s, x) (l - w {s, x) 9 {s, x) ] dxds 
Jo Jo Ju 


J J ( 2 ;) JLi (s, x) ^>1 X, 0 (s, x) , u (s, x)) dxds 
+ yyV^Ti {x)K 2 {s,x) <I >2 ^s,x,0 (s,x)) dxds 
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Let (M) denote the Hilbert space of real-valued sequences {snIneN such that (s, s) = X^neN Sn<oo. We 
set 

and 

Sy = {w € ([/;M); Vx G U,w{x) G [0,Umax]} 

Let consider the linear operator Q : D (Q) —>■ (M) x L^([/;M) where D (G) C M_|_ x (R) x L^(?7;R) 

defined by 

= \n{t)yn- / a{t,x)w{t,x)ipnix){y,ip{x))dx- / Ki{t,x)(pn{x)5{{y,(p{x))){y,ip{x))dx 
JU Ju 


G^ {t,y,z) = 0 

G is the infinitesimal generator of the evolution system E defined such as Vn G N, < s < t < ti, 

(E' (». t) {y, *-))„ = I exp (A., (r) - a (r, K, (r. .)) dr) ^ (r)) ^ (r) dr 

{s,t) iy,z) = Z 

Let also consider the operators E : D (E) PI (R+ x Sg x Sy) (R) x ([/;R) where D {E) C R_|_ x 

(R) X ([/;R) defined by 


Ei {t, y,z) = / (fn (x) (a (t, x) + K 2 (t, x) 7 {t, X, e {t, x) ,v{t,x),p {t, x)) (1- p {t, x))) dx 
Ju 

, f U AW/ / A\A (l + e)(^-^)</^n(a;) J 

-h / Ki{t,x)5 {{y,(p{x))) - dx 

Ju ^ 

I [ TX \xn r ^\-,Wn{x){y,‘p{x)) 

+ / Ki{t,x)5 {{y,(p{x))) - dx 

Ju ^ 


- / (pn{x)K 2 it,x)'y{t,x,y,v{t,x),p{t,x)){l-p{t,x))dx 

Ju 


{t, y,z)= 13 {t, x, {y, yp (x)}) ^3 (t, x, {y, ip {x )), z) 

The function E is measurable with respect to the time t. We can show as in the proof of proposition I, I. II 
that E is Lipschitz continuous with respect to (y, z) and there is a global and unique solution for the system 
dal - (135]). Note that V (t, x) G R+ x U the function f3 {t,x, .) is Lipschitz continuous and bounded on the 
set [ 0 , 1 ] using the assumption ()2.13p . ■ 

Theorem 4.2 Let consider the system ([HS]) — (|35|) and assume that T(0, .) = v (0,.). 

(i) If the functions Ki and K 2 are identically null almost everywhere on U and 

inf {a (t, x) ; t > 0, X G ?7} > 0 (36) 

then 6 is a globally exponentially asymptotically stable observer for 6. 

(ii) If the function Ki is identically null almost everywhere on U and K 2 is not identically null at least on a 
subset ofU with positive measure then 6 is a a locally asymptotically stable observers for 9. Moreover, 
if there is a positive function G (R^ x f7;R+) such that 

\l{t,x,yi,y 2 ,y 3 ) -j{t,x,zi,Z 2 ,Z 3 )\ > (t,x) \\y - z\\ (37) 

then 9 is a a globally exponentially asymptotically stable observers for 9. 
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{Hi) If the function K 2 is identically null and 

V + {1 + £ — 6) dev 


inf 


Kid ( 6 * 1 ) + a{t,x)w {t, x ); {t, x) € M+ x > > 0 . 


(38) 


then 9 is a a locally exponentially asymptotically stable observers for 9. 
(iv) If the functions Ki and K 2 are not identically null, and 
V + {l+£ - 9) dgv 


inf < Ki5 ( 9 


+ K 2^2 (t, X, 6 *^ + a {t, x) w {t, x ); {t, x) E M+ x C/| >0. (39) 

then 9 is a a locally exponentially asymptotically stable observers for 9. Moreover, if the condition dSZ]) 
is satisfied and 


K 2 {t, x) $2 9 {t, x)j > Ki {t, x) \^t, 9 {t, x), V {t, x)j , V {t, x) E x U (40) 

then 9 is a a globally exponentially asymptotically stable observers for 9. 

Proof. Let e = 9 — 9 he the error of estimation. 

(i) The error e satisfies 

dt ||e (t, .)||^ 2 ([;) /2 = - J a{t,x)w{t,x)e^ {t,x)dx 


(41) 


<-\\e{t,.)\\l^u) inf {a(s,y)} 

^ > (s,y}eR+xU 


and 


|e(t,.)||^ 2 (^) < exp I|e (0,.)||i 2 (f/) . 

We deduce the results. 

(ii) The error e satisfies 

df ||e (t, .)||^ 2 ([ 7 )/2 = — y a {t,x) w {t,x)e^dx — J {A{t,x)'Ve,'Ve) dx 
— J K 2 {t, x) ^2 (ft, X, 9 'j e^ (t, x) dx 

[ —^— —Ki{t,x)d(9^e^dx— [ a{t,x)w{t,x)e^dx 
Ju X ^ ^ Ju 


(42) 


(43) 

(44) 


< - 


j^K 2 {t,x) (1 - p) (fy{t,x,9,v,p) - 7 (t,x,9,v,p^'^e{t,x)dx 


and 


||e(L-)llL2(c/) - ■ 

Moreover, if the condition (I37p is satisfied then 

dt \\e {t ,.) ||^ 2 (^) < — y 2 q: {t, x) w {t, x) {t, x) dx — J 2 K 2 {t, x) {1 — p {t, x)) (t) {t, x) dx 


< -2 ||e(L .)||^ 2 (m inf {a {s,y)w {s,y) + K 2 {s,y) {1 - p {s,y)) {s,y)} 

^ ’ (s,y)eK+x!7 


and 


||e(t,.) 


\l^(u) 


< exp —2t 


t inf {a{s,y)w{s,y) + K 2 {s,y){l- p{s,y))C^{s,y)}] ||e (0, .)||^ 2 /m • 
{s,y)eR+xU J ^ 


Since the function C^y is positive the results holds. 
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(in) The error e satisfies 

dt \\e {t ,.) ||^ 2 ([ 7 ) /2 = — J a{t,x)w {t, x) {t, x) dx — J {^A (t, x) Ve^ (t, x) , Ve^ (t, x)) dx 
— J Ki (t, x) X, 01 {t, x ), vi {t, x)^ {t, x) dx 


< - 


V + {1 + e- 9) dev 


Ki {t, x) 5 lOi) (t, x) dx 


L 

/ a{t,x)w{t,x)e^{t,x)dx 

Ju 

r Ki {t,x) 6 ( 0 ) jd {t,0) (77Xmax {I + S - 9) - v) 


lu 


ad^-^max (1 - 0^) 

a (t) w (t) + Ki {t) 6 ( 0 )) (t, x) dx. 


-e^ (t, x) dx 


Using the condition ([551) we deduce that d is a locally exponentially asymptotically stable observers 
for 9. 


(iv) Using the same arguments developped in (ii) and {iii) we get that d is a locally exponentially asymp¬ 
totically stable observers for d. With the conditions (1371) and (I40p we are sure that even if the term 
with <l>i introduces an unstability, the stability is maintained by the term with < 1 > 2 . 


4.2 Numerical evaluation of the spatial model’s observers 

In this section we present numerical simulations in order to check the effectiveness of the observers we 
studied theoretically in the subsection 14.11 The duration of the simulations is the year. Parameters are 
taken following m and generalize those in the subsection 13.21 We assume a anisotropic radial control 
strategy u given for every (f, x) € M_|_ x U, 

u{t,x) = sin^ (^\\M (x - xo)|P^ sin^ {t - exp (^-uj 2 {t - • (45) 

The matrix M introduces a geometrical anisotropy in the distribution of u. The functions a, 13 and 7 are 
taken such as 

a {t, x) = Pi {t, x) -|- biqi (x) (1 — cos (cit)) {t — di)^ , (t, x) G M+ x U (46) 

f3 {t, X, y) = 62^2 (a:) (1 - cos (c 2 t)) {t - d 2 f P 2 {y), (t, x, y) G M+ x U x M (47) 

and 

l{t,x,yi,y 2 ,y 3 ) = dads (x) (1 - cos (cat)) {t - d^f {yi - ^^ 3 ) 2 / 2 , {t,x,y) G M+ x U x M^. (48) 

Pi is a nonnegative function of the time, p2 is a positive function and for every i in {1,2,3}, qi is a real 

anisotropic radial nonnegative function of the space variable as the control u. For every i in {1,2,3}, 
bi, Ci, and dj are positive coefficients corresponding respectively to the maximal amplitude, the pulsation 
and the global maximun of a, (3 and 7 . For every i in {1,2,3}, qi is a [0, l]-valued function of the space 
variable. The matrix A in equation ([T]) is assumed to be equal to 10“^/ because the diffusion is relatively 
slow. The functions Ki and K 2 are taken constant and are similar with ki and /c 2 defined in subsection 
[321 The rest of the parameters are the same described in the table [T1 The initial conditions are also taken 
following the subsection 13.21 and are constant with respect to the space variable in order to easily fulfill the 
boundary conditions (l2|). The following table specifies some assumed values of additional parameters: For 
every i G {1; 2; 3}, 
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Parameters 

Values 

Source 

Parameters 

Values 

Source 

Qi {x) 

Xi 

sin^(||Mi(a;-Xi)||^)+l 

2 

0 

Assumed 

Assumed 

M,M, 

randomly generated 

5 * rand (3,3) 

Assumed 


Table 2: Simulation parameters for the spatial model 


In the remainder of the subsection we present simulations of the observers and their estimation absolute 
relative error dynamics. To keep the representation simple we just plot the spatial minimum values, the 
spatial average values and the spatial maximum values. Instead of the estimation relative error we consider 
the estimation absolute relative error to avoid compensations when averaging. The initial condition p{0 ,.) 
is always taken equal to 0(0,.). Indeed, following simulations in subsection 13.21 the performances of the 
observers decrease with respect to p (0,.) and on the other hand, it seems realistic to take p (0,.) less than 
0 ( 0 ,.). 

4.2.1 Simulations with 0 (0 ,.) = v (0,.) = 0.05 

The figures [9] and [10] show the behaviour of the observers when the observation is started early. As we can 
see the global performance of each observer is relatively good but the case ki = 0 and k 2 = 10^ seems to be 
the better one. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and k 2 = 0 


ki = 10^ and k 2 = 10^ 


Figure 9: Inhibition rate and observers for 0 (0 ,.) = v (0,.) = 0.05. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 ki = 10^ and /c 2 = 10^ 

Figure 10: Relative absolute estimation error for 9 (0 ,.) = v (0,.) = 0.05. 
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4.2.2 Simulations with 6 (0,.) = 0.05 and v (0,.) = 0.5 

The figures fTTl and fT^ correspond to an observation starting at the beginning of the disease but with a well 
developped berry. We have better performances when k 2 = 10^. With ki = 10^ and ^2 = 0 we notice an 
unstability behaviour in figure [T2] although the error takes the null value. That may be due to the constant 
sign of the function 4>i. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and ^2 = 0 ki = 10^ and k 2 = 10^ 

Figure 11: Inhibition rate and observers for 0 (0,.) = 0.05 and v (0,.) = 0.5. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 ki = 10^ and /c 2 = 10^ 

Figure 12: Relative absolute estimation error for 9 (0,.) = 0.05 and v (0,.) = 0.5. 
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4.2.3 Simulations with 6 (0,.) = 0.75 and v (0,.) = 0.05 

In the figures [TJ] and [T4] the inhibition rate is already high while the berry is little. As expected the worst 
estimation corresponds to ki = k 2 = 0. The best estimations are still made when k 2 = 10^. However, when 
ki = 10^, the same unstability behaviour in figure [T2] appears. 


ki = k2 = 0 


ki = 0 and k 2 = 10^ 


ki = 10^ and /c 2 = 0 ki = 10^ and /c 2 = 10^ 

Figure 13: Inhibition rate and observers for 9 (0,.) = 0.75 and v (0,.) = 0.05. 
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fci = /C2 = 0 


ki = Q and k 2 = 10^ 


ki = 10^ and /c 2 = 0 


ki = 10^ and /c 2 = 10^ 


Figure 14: Relative absolute estimation error for 9 (0,.) = 0.75 and v (0,.) = 0.05. 
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4.2.4 Simulations with 6 (0,.) = 0.75 and v (0,.) = 0.5 

For the figures [TJ] and [TU] the inhibition rate is already high and the fruit is relatively mature when the 
observation starts. The worst case are still when ^2 = 0 and the best are with k 2 = 10^. 


= ^2 = 0 


fei = 0 and k 2 = 10^ 


ki = 10^ and A :2 = 0 ki = 10^ and /c 2 = 10^ 

Figure 15: Inhibition rate and observers for 9 (0,.) = 0.75 and v (0,.) = 0.5. 
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fci = /C2 = 0 


ki = d and k 2 = 10^ 


ki = 10^ and /c 2 = 0 


ki = 10^ and /c 2 = 10^ 


Figure 16: Relative absolute estimation error for 6 (0,.) = 0.75 and v (0,.) = 0.5. 
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5 Conclusion and prospects 


We have studied four observers for the dynamics of anthracnose disease. Our constructions were based on 
the models given in [12]. We consider two models: a spatial version and a non spatial one stated respectively 
by the equations ([I|) — (l6|) and ([7|) — (jlOp . Some changes have been made in the original form of those models 
in order to be more realistic. Precisely, we have changed the equations modelling the dynamics of the berry 
volume (v) and the rot volume (vr)- Indeed, before the launch of the disease, the berry has to reach a 
minimal volume. On the other hand the rot volume shall never be greater than total volume. We found that 
the system describing the disease dynamics can be viewed as the natural observer of the system even the 
performances are poor. We have studied theoretically two observers which gives fairly good results looking 
at the relative error. Several simulations have been carried out to assess the effectiveness of our proposed 
observers. That work has been done for the spatial and the non spatial models. 

Although proposed observers seem to display fairly good results, we think that a stochastic approach 
could be introduced in the estimation process. Indeed, the parameters of the models are submitted to several 
stochastic variations. Moreover, the measures of v and Vr also contains errors. A mean to take into account 
those errors could consist in adding noises in our models and proceed to the stochastic filtering by different 
existing tools. That aspect is studied in on going works of some of the authors of the current paper. 
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